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For a random walk on a network, the mean first-passage time from a node i to another node j chosen 
stochastically according to the equilibrium distribution of Markov chain representing the random walk is 
called Kemeny constant, which is closely related to the navigability on the network. Thus, the configuration 
of a network that provides optimal or suboptimal navigation efficiency is a question of interest. It has been 
proved that complete graphs have the exact minimum Kemeny constant over all graphs. In this paper, 
by using another method we first prove that complete graphs are the optimal networks with a minimum 
Kemeny constant, which grows linearly with the network size. Then, we study the Kemeny constant of a 
class of sparse networks that exhibit remarkable scale-free and fractal features as observed in many real-life 
networks, which cannot be described by complete graphs. To this end, we determine the closed-form solutions 
to all eigenvalues and their degeneracies of the networks. Employing these eigenvalues, we derive the exact 
solution to the Kemeny constant, which also behaves linearly with the network size for some particular cases 
of networks. We further use the eigenvalue spectra to determine the number of spanning trees in the networks 
under consideration, which is in complete agreement with previously reported results. Our work demonstrates 
that scale-free and fractal properties are favorable for efficient navigation, which could be considered when 
designing networks with high navigation efficiency. 
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Kemeny constant, defined as the mean first- 
passage time from a node i to a node j selected 
randomly according to the equilibrium distribu- 
tion of Markov chain, is a fundamental quantity 
for a random walk, since it is a useful indicator 
characterizing the efficiency of navigation on net- 
works. Here we prove that among all networks, 
complete graph is the optimal network having 
the least Kemeny constant, which is consistent 
with the previous result obtained by a different 
method. Then, we study the Kemeny constant of 
a class of sparse scale-free fractal networks, which 
are ubiquitous in real-life systems. By using the 
renormalization group technique, we derive the 
explicit formulas for all eigenvalues and their mul- 
tiplicities of the networks, based on which we de- 
termine the closed-form solution to the Kemeny 
constant. We show that for some particular cases 
of networks the leading scaling of Kemeny con- 
stant displays the same behavior as that of com- 
plete graph, indicating that for networks with 
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many nodes, navigation performance similar to 
that of complete graph can be obtained by net- 
works with many few links. Finally, to show the 
validity of the eigenvalues and their degenera- 
cies, we also use them to count spanning trees 
in the networks being studied, and recover previ- 
ously reported results. Our work is helpful for the 
structure design of networks where navigation is 
very efficient. 



I. INTRODUCTION 

In the past decade, many research reports demon- 
strated that complex networks are ubiquitously present 
in nature and social sciences^—. Today, complex net- 
works have been recognized as an emerging subject and 
useful tool for studying complex systems^. An outstand- 
ing issue in network science is to understand the behav- 
ior of various dynamical and stochastic processes taking 
place on different networks with rich topologies, which is 
not only of theoretical interest but also of practical rele- 
vance^. Among diverse dynamics, random walka&"— are 
a paradigmatic process and have been extensively stud- 
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ied 9 - — due to their wide range of applications in science 
and engineering^—. 

A fundamental quantity for random walks is mean 
first-passage time (MFPT)22. The MFPT from a source 
node i to a destination node j, denoted by Fij, is the ex- 
pected time for a walker starting from i to reach j for the 
first time. Many other quantities related to random walks 
can be expressed in terms of, or encoded in, MFPT. For 
example, for the trapping problem^ — a particular case 
of random walks with a trap located at a given node ab- 
sorbing walkers visiting it — the trapping efficiency quan- 
tified by average trapping time is defined as the mean 
of node-to-trap MFPT over source nodes. As a quanti- 
tative indicator of trapping efficiency, average trapping 
time has been studied on various networks with distinct 
topological properties, with an attempt to uncover the 
impact of network topology on the trapping process. Ex- 
amples include classic fractals^—, non-fractal scale-free 
networks 3 ^—, fractal scale-free networks^—, modular 
scale-free network o 39 ' 40 , and so on. 

In addition to the trapping problem, another interest- 
ing quantity derived from MFPT is the Kemeny constant, 
which is defined as the expected time for a walker start- 
ing from a node i to another node j selected randomly 
from all nodes according to the stationary distribution 
of Markov chain representing the random walk 41 -. Since 
the Kemeny constant is equal to the sum of reciprocals of 
one minus each eigenvalue of the transition matrix other 
than eigenvalue one 4 ^, it is also called eigentime iden- 
tity^ 3 -, which is independent of the starting point. It can 
be explained in terms of navigability^ 4 -, measuring the 
efficiency of navigation. Thus, it is of theoretical and 
practical significance to study the Kemeny constant of 
complex networks. In contrast to the average trapping 
time, little attention has been paid to the Kemeny con- 
stant, and many questions about this relevant quantity 
are still open. 

In this paper, we address a specific question: of all 
networks, which have a minimum or almost minimum 
Kemeny constant? It is of theoretical relevance and prac- 
tical importance, and is particularly useful for network 
design. By making use of the connection between the 
Kemeny constant and eigenvalues of the transition ma- 
trix, we first show that among all undirected networks 
complete graphs have the minimum Kemeny constant 
growing linearly with the network size, which agrees with 
result reported before 4 ^. However, since most real net- 
works are spare in the sense that their average degrees 
are much less than that of complete graphs. Moreover, it 
was shown that many real networks exhibit the striking 
scale-free 4 ^ and fractal 4 ^ properties. Here we study the 
Kemeny constant for a class of sparse scale-free fractal 
networks 4 ^, denoted by H(q,n), which are parameter- 
ized by a positive integer q > 2. Using the decimation 
technique, we derive explicit expressions for all eigenval- 
ues and their multiplicities of the transition matrix of 
-ff (<7, n), based on which we derive a closed- form formula 
for the Kemeny constant. The obtained result shows 



that, for q > 3, the Kemeny constant also scales lin- 
early with the network size, displaying the same scaling 
as that of complete graphs. Finally, we use the eigen- 
values of H(q,n) to determine the number of spanning 
trees and recover previous results, which also validates 
our computation of eigenvalues for such networks. 

II. FORMULATION OF THE PROBLEM 

Consider discrete-time unbiased random walks^ on a 
connected undirected network G with N nodes and E 
edges, where the N nodes are labeled by 1, 2, 3, . . . , N, 
respectively. The connectivity of G is encoded in its ad- 
jacency matrix A, whose entry a.^ = I (or 0) if nodes i 
and j are (not) connected by an edge. The degree of node 
i is defined to be di = Y^j=i a iji an d the diagonal degree 
matrix of G, denoted by £), is defined as: the ith diagonal 
element is di, while all non-diagonal elements are zero. 
Thus, the total degree of all nodes is K = 2E = J2jLi di, 
and the average degree is (d) — 2E/N. 

For the unbiased random walks considered here, at ev- 
ery time step the walker starting from its current lo- 
cation jumps to each of its neighbors with an identi- 
cal probability. Such a stochastic process is character- 
ized by the transition matrix T = D~ 1 A, whose entry 
Uj — o-ij/di presents the probability of moving from i 
to j in one step. Note that the transition matrix is also 
called the Markov matrix, since a random walk is an er- 
godic Markov chai n 41 i 43 , whose stationary distribution 
7r = (tti, 7T2, . . . , 7r/v) T is a unique probabilistic vector 

satisfying 7r$ = di/ K, Yli=i n i — 1; an d k t T = tt t . 

In general, the transition matrix T is asymmetric ex- 
cept for regular networks. So, we introduce a matrix P 
similar to T by 

P = D'iAD-i = D^TD-i, (1) 

where is defined as follows 4 ^: the ith diagonal entry 
is l/\/di, while all non-diagonal entries are equal to zero. 
It is evident that P is real and symmetrical and thus has 
the same eigenvalues as T. 

Various interesting quantities related to random walks 
are encoded in matrix P, and thus can be deduced from 
P. For example, the MFPT F^ for a walker starting from 
a node i to another node j can be expressed in terms of 
the eigenvalues and their orthonormalized eigenvectors of 
matrix P. Let Ai, A2, A3, . . ., Xn be the N (real) eigen- 
values of matrix P, rearranged as 1 = Ai > A2 > A3 > 
• ■ ■ > Ajy > — 1. Obviously, Ai + A 2 + A 3 + • • • + X N = 0. 
And let ipi, -02, V>3j ■ • ■> V'JV denote the corresponding 
normalized, real-valued and mutually orthogonal eigen- 
vectors, where fa = (fai,fa2, ■ ■ ■ , AnV ■ Then, 

which has been derived by various methods, e.g., gener- 
ating functions 8 and spectral graph theory 4 ^. 



Another interesting quantity for random walks is the 
Kemeny constant^ 1 -, also referred to as eigentime iden- 
tity^, defined by 
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Equation (j4]) indicates that F is the average of MFPT 
Fij from node i to node j randomly chosen from all 
nodes accordingly to the stationary distribution. Since 
F is independent of the starting point i, it is thus called 
the Kemeny constant. This quantity is closely related 
to user navigation through World Wide Web and can 
be explained as the mean number of edges the random 
surfer needs to follow before reaching his/her final desti- 
nation^. 

Note that the Kemeny constant is a global spectral 
characteristic of a network ^ 41 ' 43 . It has been show n 8 ' 41 ' 43 
that 
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which holds for any connected undirected network. Ke- 
meny constant can be used as a measure of navigation 
efficiency 4 ^ and is thus relevant to Web design: small F 
can save user's effort by directly guiding them to most 
desirable Web pages. Since F depends on the network 
structure, a question arises naturally: what is the opti- 
mal network structure that has a minimum Kemeny con- 
stant? Intuitively, the globally optimal network should 
be fully connected. In the sequel, we will address this 
question analytically by using the spectra of graphs and 
method of Lagrange multipliers, and show that the eigen- 
value treatment arrives at the intuitive and natural con- 
clusion. 



III. GLOBALLY OPTIMAL NETWORK 

We now solve the aforementioned problem by using 
the method of Lagrange multipliers. We consider the 
following function: 



f(x 2 , ■ ■ ■ ,x N ) 
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subject to the condition of 1 + x 2 + ■ • • + Xjv = on the 
domain 1 > x% > X3 > ■ • • > xn > —1. Applying the 
method of Lagrange multipliers to Eq. ([5]), the minimum 
of the function f{x%, ■ • • ,xn) is attained at x 2 = £3 = 
• • • = xn = —1/(N — 1). Thus, the minimum Kemeny 
constant among all networks of size N is 



(N-l) 2 
N 



(6) 



As one can see, the above set of identical eigenval- 
ues A2 = A 3 = . . . = Ajy = —1/(N — 1) corresponds to 
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FIG. 1. (Color online) Construction of the networks. One 
can obtain the next generation of the network family through 
replacing each edge of the present generation by the clusters 
on the right-hand side of the arrows. 



those of the transition matrix of a complete graph with 
A nodes, in which every node is connected to all oth- 
ers. That is, among all connected undirected networks, 
complete network is optimal in the sense that it has the 
absolute minimum Kemeny constant, which can not be 
achieved by other networks. We note Eq. © has been 
previously proved by using a different approach 43 . It is 
obvious from Eq. © that for large N, F m i n scales linearly 
with N. 

Although a complete network has the minimum Ke- 
meny constant, it is dense, not like most real networks. 
Extensive empirical research has shown 1 -^ that most real 
networks are sparse, which have a small average de- 
gree and display simultaneously scale- free^ and fractal 4 ^ 
properties. Below, we will study the Kemeny constant for 
a family of fractal scale-free sparse networks 4 ^ parame- 
terized by a positive integer q > 2. We will show that 
for those networks corresponding to q > 3, the dominat- 
ing term of their Kemeny constant grows linearly with 
the network size A, which is similar to the linear behav- 
ior of complete graphs but with a different prefactor 2^1 
larger than 1. Thus, when q > 3, the networks are called 
"suboptimal networks" . 



IV. SUBOPTIMAL SCALE-FREE FRACTAL 
NETWORKS 



In the section, we study random walks on a class of 
scale-free (heterogenous) fractal networks 4 ^, the prop- 
erties of which are dominated by an integer parameter 
q > 2. 
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These resultant networks display some remarkable fea- 
tures observed in various real-life systems. They have 
a power-law degree distribution with exponent 7 = 2 + 
ln2/lng, and have a fractal dimension /b = ln(2g)/ln2, 
both of which show that they are fractal scale-free net- 
work o 47 ' 49 i 50 . Moreover, they are "large- world" , with 
both of their diameter and average distance increasing 
as a power function of the network siz o 47 i 49 ' 50 . We note 
that this "large- world" phenomenon was also observed in 
the global network of avian influenza outbreak o 52 ' 53 . 



B. Spectrum of the transition matrix 

After introducing the construction method and prop- 
erties of the networks, we proceed to determine the eigen- 
values and their degeneracies of the transition matrix for 
the network family, which will be useful for evaluating 
the Kemeny constant. 



1. Eigenvalue spectrum 



FIG. 2. (Color online) Illustration for the network H(2,5). 



A. Construction and properties 

The scale-free fractal networks studied in~ are con- 
structed in an iterative manner as depicted by Fig. Q] 
Let H(q,n) (q > 2 and n > 0) denote the networks af- 
ter n iterations. Initially (n = 0), H(q,0) is an edge 
connecting two nodes. For n > 1, H(q,n) is gener- 
ated from H(q, n — 1) by replacing every existing edge in 
H(q,n—1) with the clusters on the right-hand side of the 
arrows in Fig. [TJ When q = 2, H(q, n) is reduced to the 
(2, 2)— flower, which is a special case of the (u, v)— flowers 
(u > 1, v > 2) presented i n 49 ' 50 . Figure [5] illustrates the 
network H(2, 5) for the particular case of q = 2, which 
is called diamond lattice and was introduced more than 
thirty years agc£k 

The iterative construction allows to precisely analyze 
relevant properties of the networks. At each generation 
n-i {tii > 1), the number of newly introduced nodes is 
Kij = q(2q) ni ~ 1 . Then, the network size (number of 
nodes) N n of H(q, n) is 



N„ 



g(2g) n + 3g - 2 
2q- 1 



(7) 



Let E n denote the number of total edges in H(q,n). It 
is easy to verify that 



E n = qE n ^ = (2q) r 



(8) 



Let di(n) be the degree of node i in H(q,n) that was 
generated at iteration rii (rii > 1). Then, we have di(n + 
l)=qdi{n) = 2q n - n K 



Let A n and D n denote, respectively, the adjacency ma- 
trix and diagonal degree matrix of H{q 1 n). The entry 
A n (i,j) of A n is defined as: A n (i, j) = 1 if nodes i and j 
are directly connected in H(q, n), A n (i, j) = otherwise. 
Then, the transition matrix of H(q,n), denoted by T n , 
is defined as T n = (D n )~ 1 A n . Thus, the entry of T n is 
T n (i,j) = A n (i,j)/di(n). 

We now begin to compute the eigenvalues of T n . Since 
T n is asymmetric, we introduce a real symmetric matrix 
P n similar to T n by 



(D n )-iA n (D n y 



(D n )*T n (D n y 



(9) 



Thus, we reduce the problem of determining eigenvalues 
of T n to calculating those of P ni with entries P n (i,j) = 

Next, we apply the decimation method^— to com- 
pute the eigenvalues and their multiplicities of P n . Let 
us consider the eigenvalues of matrix P n +i- For this pur- 
pose, we use a to denote the set of nodes belonging to 
H(q,n), and (3 the set of nodes created at iteration n + 1. 
By definition, P n +i has the following block form: 



Pn+l — 



' p 






P a ^~ 








Pf},a 



(10) 



where we have used the fact that P a . a and Pp.p are two 
zero matrices with orders being N n x N n and (N n+ i — 
N n ) x (N n +i — N n ), respectively. The reason for P QjC t and 
Pp,p being zero matrices lies in the iterations and inser- 
tion constraint in the construction process of networks, 
as shown in Fig. [1] 

Suppose Aj(n + 1) is an eigenvalue of P n +i, and u = 
(u a ,up) T is its associated eigenvector, where u a and up 
correspond to nodes in a and /3, respectively. Then, the 
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eigenvalue equation for matrix P„+i has the following 
block form: 



P a ,p ' 






Pft,a 







= \i(n + l) 

which can be rewritten as 

P a ,pup = \ i (n + l)u a , 



Up 



(11) 



(12) 



Pp, a u a = \i{n + l)up. (13) 
In the case of \{n + 1) ^ 0, Eq. (IT3")) gives rise to 



Xdn + 1 



(14) 



Inserting Eq. flU) into Eq. ([12]) yields 
1 



A 



I _u i\ P <*,P p P,»u a = A t (n + 1)m q . (15) 



Thus, the problem of determining the eigenvalues for ma- 
trix P n+ i of order -/V, i+1 x iV n+1 is reduced to calculating 
the eigenvalues of matrix P a ,pPp, a with a smaller order 



In Appendix lAl we prove that 



1 1 



(16) 



where /„ is the identity matrix of order N n x N n , same 
as that of P n . Thus, we have related P a ,pPp, a to P n , 
which allows to express the eigenvalues of matrix P n +\ 
in terms of those of P n . 

Plugging Eq. (fTo) into Eq. (IT51 gives 



J n + \p n ) u a = Ai(n + l)u„, (17) 



Ai(n + 1) V2 2 
that is, 

P„u Q = {2[A 4 (n + l)] 2 -l} 



(18) 



Hence, if Aj(n) is an eigenvalue of P„ corresponding to 
the eigenvector u a , then Eq. (ITS)) implies 



X t (n) = 2[X t {n + l)} 2 - 1. 



(19) 



Solving the quadratic equation (|T9|) in variable Aj(n + 1) 
gives two roots: 



x , / Aj(n) + 1 / Ai(n) + 1 
Ai,i(n+l) = V ^ ' A *,2(«+l) = -y ^ , 

which relate Ai(n+1) to \i(n), with each Ai(n) producing 
two eigenvalues of P n +i. 



2. Multiplicities of eigenvalues 

To determine the multiplicities of eigenvalues, we first 
calculate numerically the eigenvalues for those networks 
of small sizes. For n = 0, the eigenvalues of Po are 1 and 
— 1 ; while for n = 1 , the eigenvalues of matrix Pi corre- 
sponding to H q i are 1,-1, and with multiplicity q. For 
n > 2, we find that the eigenvalue spectrum displays the 
following properties: (i) Eigenvalues 1 and —1 are present 
at all generations, both having a single degeneracy, (ii) 
All eigenvalues of generation rii will always exist at its 
subsequent generation m + 1, and all new eigenvalues at 
generation m + 1 are just those generated via Eq. (f2T)]) by 
substituting Aj(n) with A,(nj) that are newly added to 
generation nf, moreover, each new eigenvalue keeps the 
degeneracy of its father, (iii) Except for 1 and — 1, all 
other eigenvalues are derived from generated 0. Thus, all 
that is left to determine are the degeneracy of and the 
multiplicities of its offsprings, based on property (ii). 

Let D™ ul (A) denote the multiplicity of eigenvalue A of 
matrix P n , We now find the number of eigenvalue of 
P n . Let r(M) denote the rank of matrix M. Then, the 
degeneracy of eigenvalue of P„+i is 



^+i(A = 0)=7V„ +1 -r(P„ +1 ). (21) 

In Appendix [B] we derive the exact expression for the 
degeneracy of eigenvalue of P„ as 



D™\\ = 0) = 



0, 

(?~l)(2g)" 
2q-l 



n = 0, 
n > 1. 



(22) 



Since every eigenvalue in P n keeps the degeneracy of its 
father, the multiplicity of each first-generation descen- 
dants of eigenvalue is ^"^-i — ~> tne multiplicity 
of each second-generation descendants of eigenvalue is 

(g ~ 1} 2g-i" 2+g ' and 80 on - Thus ' for Pn ( n ^ X )' the total 
number of eigenvalue and all of its descendants is 



K cod (o) = £ 



(q-l)(2q¥+q 



i=i 

q(2q) n - q 



2q-l 



x 2 r " 



2,-1 < 23 > 

Summing up the multiplicities of eigenvalues obtained 
above, we have 



q{2q) n + 3g - 2 
2q- 1 



= N n . 
(24) 

Thus, we have found all the eigenvalues of P n . 

Thus far, we have given a recursive expression for 
eigenvalues of the transition matrix for H(q, n), as well as 
the multiplicity of each eigenvalue. In fact, for the special 
network family H(q, n) considered here, their eigenval- 
ues can be expressed in an explicit way, as shown in Ap- 
pendix [C] These exact expressions contain more concrete 
information about eigenvalues than that of the recursive 
relation provided by Eq. (|2T))) . 
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C. The Kemeny constant 

We proceed to use the obtained eigenvalues and their 
multiplicities to determine the Kemeny constant for the 
fractal scale-free networks H(q, n). 

Let Fij(n) denote the MFPT from node i to node j 
in H(q, n). Let it = (7r 1; 7r 2 , . . . , tt./v„) T represent the sta- 
tionary distribution 42,43 for random walks on H (q, n). It 
is easy to verify that iri = di(n)/(2E n ). Then, the Ke- 
meny constant for H(q,n), denoted by F n , is given by 



which indicates that 

E ^T- 4 E 



1 



1 



(31) 



Thus, we have 



N„ 



i 



Ai(n) 



(25) 



where we have assumed that Xi(n) = 1. 

Let L n be the normalized Laplacian matri x 48 ' 58 of 
H(q, n), defined as L n = 7„-P„ = /„- (D n )?T n (D n )~? . 
Let (Ti, (72, <t 3 , (jjv„ be the iV n real non- negative 
eigenvalues of matrix L n , rearranged as = a± < oi < 
cr 3 < ... < (jjv„ = 2. Note that for any i, <Ti{n) 
1 — Aj(n). Then, 



E — 

' a An) 



(26) 



Next, we explicitly evaluate this sum. 

Let n be the set of all the N n — 1 nonzero eigenval- 
ues of matrix L n , i.e., fl n = {a2(n), 03(71), . . . , crN n (n)}, 
in which the distinction of the elements has been ig- 
nored. It is clear that Sl n (n > 1) includes 1, 2, and 
other eigenvalues generated by 1. Then, Sl n can be clas- 
sified into three subsets represented by Sln\ fin'' and 

(3) 

Sl n , respectively. That is, Sl r , 



where Sln^ consists of eigenvalue 1 with multiplicity (2q) n = ( 2< 3~ 1 ) jv "~ 3l ?+ 2 anc i 4™ = 



Sin U Sin U Sin • 



F „ = 4F„_ 1+ *i^»l±j -| (32) 

Considering the initial condition Fq = \ , Eq. (|32p can be 
solved by distinguishing two cases: q — 2 and 9 > 3. For 
9 = 2, ' 

F n = -L (6n x 4" + 4" +1 + 5) ; (33) 
lo 

while for q > 3, 

_ 3(q - l){2q) n+1 - 2q(2q - 1)4" + Aq 2 - 1L? + 6 

U_ 6(g-2)(2g-l) ' 

(34) 

We now express i 7 ^ as a function of the network size 
N n . For the case of q — 2, we have 4 71 = |JV„ — 2 and 
ri = log 4 (|iV n — 2), which can be easily obtained from 
Eq. ([?]). Thus, Eq. (1331) can be rewritten in terms of N n 
as 

F n = QiV n - I J log 4 f |jV„ - 2) + 67Y„ - 3 . (35) 
Similarly to the case of q = 2, for q > 3 we can derive that 



(g-l)(2g)"+g (2) 



2,-i > n » contains only eigenvalue 2 with single Inserting these re l a tions into Eq. 



(3) 

degeneracy, and f2^ includes those eigenvalues gener- 
ated by 1. Obviously, E ffi(n)e ^) ^TJ = ^"^T^ 



(2g-l)JV„-3g+2 
1 

leads to 



Pn = 



-N n 



22q 2 -41q+ 18 



and ^ , x ^„(2) — 7 — r — t;. 

For each eigenvalue (?i(n— 1) in fi„_Li, ^ generates two 

(31 

eigenvalues, <7i,i(n) and <Ti^{n), belonging to Sl n , via 
the following equation: 



3(9-2) 



6(9-2)(2g 
\2q-\)N n 



1) 

3g 



1 lo S2 



(36) 



2[cr l (n)] 2 - 4a z (n) + <j t (n - 1) = , 



(27) 



which can be easily obtained from Eq. (fl~9|) . According 
to Vieta's formulas, we have 



Equations ([33)) and (|3"6"|) show that, for large networks, 
i.e., iV ff — > co, the leading terms in the Kemeny constant 
for -ff (g, n) follow distinct scalings for the two cases of 
q = 2 and q > 3. For the former case, 



F n ~ AL InW, 



Oi,i( n ) + <r ii2 (n) 



and 



Then 



<7»,i(n) X Ci,2(n) 



cri(n - 1) 



1 



1 



Oi,x{n) o-ji(n) (Ti(n-1) ! 



(28) 



(29) 



(30) 



while for the latter case, 



F n ~ N n . 



(37) 



(38) 



It should be stressed that, for q > 3, the linear growth 
of the Kemeny constant for sparse H(q,n) is similar to 
that for dense complete graphs. In this sense, we say 
that the scale-free fractal networks H(q, n) (for q > 3) 
are suboptimal networks. 
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Finally, we note that the validity of the above-obtained 
results for the Kemeny constant depends on that of the 
eigenvalues. In Appendix |D] to verify that our computa- 
tion on the eigenvalues of the transition matrix is correct, 
we compute the number of spanning trees in H(q, n) us- 
ing the connection between spanning trees and the eigen- 
values of the normalized Laplacian matrix, which repro- 
duces previously obtained results for the number of span- 
ning trees in H (q, n)2». This indeed validates our com- 
putations. 



V. CONCLUSIONS 

To conclude, by applying a different method from be- 
fore, we have shown that among all networks, complete 
graphs are optimal in the sense that they have the min- 
imum Kemeny constant, which scales linearly with the 
network size. However, due to the limitations of cost 
and other factors, complete graphs are not typical in the 
real world, where most systems are described by sparse 
networks displaying simultaneously scale- free and fractal 
properties. To explore networks having such properties 
of real systems but with a small Kemeny constant, we 
have also studied random walks on a scale-free fractal 
network family, characterized by a positive integer pa- 
rameter. We have derived explicit formulas for all eigen- 
values and their multiplicities of the transition matrix for 
such networks. Using these results, we have further deter- 
mined the eigentime identity, which also displays a linear 
growth with the network size as complete graphs. At last, 
to test our results on the spectra of the transition matrix 
for the scale-free fractal networks, we have computed the 
number of their spanning trees, which is consistent with 
those previously obtained results by means of a different 
method. Our work may have important implications on 
the design of complex networks with efficient navigation. 
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as 



Pa,0 = 



( M i \ 
\ M nJ 



(Al) 



where each Mj is a column vector of order N n+ i — N n . 
Since P a> p = Pj a , we have P@ ja = (Mi M 2 ...M N J. 
Then, the entry R n (i,j) = Mj Mj of R n can be deter- 
mined by distinguishes two cases, as follows. 
If i = j, then the diagonal element of R n is 

R n (i,i)=MjM i 

= P n +i (i, k)P n+1 (k, i) 



E 

fce/3 



A n+ i(i,fc) 



di(n + l)dk(n + 1) 



U(n + i) ^ 2 



1 

2 ' 



(A2) 



where i ~ fc stands for that nodes i and fc are adjacent 
in H q>n+ i and the fact d k (n + 1) = 2 was used. 
If i ^ j, the non-diagonal element of Rn is 

R n (i,j) = M?M j 

= ^2Pn+l(i,k)P n+1 (k,j) 



E 



1 d k (n + l)^(n+l)<fy(ra+l) 



A n+1 (i,k) 

2^(71)^(71) 2 n[hJ) 



Uhj), (A3) 



where the relations di(n + 1) = qdi{n) and dj(n + 1) = 
qdj(n) have been used. Thus, Eq. (fT6| is proved. 
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Appendix A: Proof of equation (|16|l 

In order to prove P a ,pPp,a = \ln + \Pm it suffices to 
show that their corresponding entries are equal to each 
other. For simplicity, let Q n = \l n + \P n , the entries of 
which are: Q„(i,i) = \ for % = j and Q n {i,j) — hP n (i,j) 
otherwise. Let R n — P a ,0Pp,a, whose entries R n (i,j) can 
be determined as follows. Note that P a ,p can be written 



Appendix B: Derivation of equation {[22) 

Equation ([2~Tj) shows that in order to determine 
D™\(\ = 0), we can alternatively compute r(P n+ i). 
Obviously, r(P n+ i) = r(P a , p ) + r{Pp, a ) = 2r(P /3)Q ), 
where r(P a ^) — r(Pp jCl ) is used. 

We next determine r(Pp. a ). As defined in Ap- 
pendix El Pp ia = (Mi M 2 ■••M Nn ), where M l = 
(M M , M 2 ,i, ■■■ , M Nn+1 - Nn ,i) T . We first show that Mi, 
M2, ■ ■ ., and Mjv„ are linearly dependent. To this end, 
we label one node in H(q,0) by 1. Then, all nodes in 
H(q, n + 1) can be categorized into 2" +1 + 1 classes 
Cj (i = 0, 1, . . . , 2 n+1 ), according to the shortest dis- 
tance between any of them and node 1. Let d(x,y) 
denote the shortest distance from node x to node y 
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in H(q,n + 1), and Cj be the set of nodes denned by 
Ci = {u\d(l,u) = i}. Then, the set a of old nodes be- 

2" 

longing to H(q,n + 1) is a — C^i, and set j3 of new 

i=0 

2"-l 

nodes born at generation n + 1 is (3 = C^-i- It is 
not difficult to verify that 



E [i-^ 1 ' 



o. 



(Bl) 



Therefore, column vectors Mi, M2, . . ., and Mjv„ are 
linearly dependent. 

We next demonstrate that the N n — 1 vectors M2 , M3 , 
. . ., and are linearly independent. Let 



v = (v 2 ,v 3 , . . . ,v Nrl+1 _ Nrl ) T 



E 

i=2 



(B2) 



Suppose that v = 0, and we will prove that for an ar- 
bitrary ki, ki = always holds. Note that for a node 
w G P, it has two neighbors, denoted by w pre and w suc : 
w pic is the neighbor near by node 1, while w suc is the 
other neighbor farther from node 1. If w G C\, then 
v w = k w + k Wsuc . Since, in this case, k Wplo is just 
node 1, v w = leads to k Wsuc = 0. Let's consider a node 
j € C2 • By construction, there exists a node w £ Ci such 
as w suc = j- Thus, if Vj = 0, kj = holds for j e C2. 
In a similar way, we can prove that for any node i £ a 
(i 7^ 1), if v = (0,0,..., 0), then fc, = 0. In this way, 
we have proved that Mi, M2, . . ., and Mjv„ are linearly 
independent. Therefore, r(P (3 , Q ) = N n - 1 = g(2g 2 ) g "^~ 1 . 
Inserting this formula into Eq. (I21[) . we can obtain the 
degeneracy of eigenvalue of P n , given by Eq. (|22|) . 



Appendix C: Explicit expression of eigenvalues of P n 

Notice that all the eigenvalues of P n lie between —1 
and 1. Let A n be an eigenvalue of P n (n> 0). Then, one 
can assume that A n = cos(# rl ) with < 6 n < n. Thus, 
from Eq. (|19[) . we have 



cos 2 (d n ) = cos 2 (0 n _i/2), (CI) 

which implies that 

cos(0 n ) = ±cos(d„_i/2). (C2) 

Under the constraint of < 9 n < ir, Eq. (|C2[) means 

0„ = n _ x /2 or 9 n = 7r-e n _i/2, (C3) 

which enable us to express the eigenvalues for P n (n > 0) 
in an explicit way as follows. 



We use an ordered pair (A„,_D™ ul ) to represent an 
eigenvalue A„ = cos n of matrix P n and its correspond- 
ing multiplicity D™ ul . Let $ n denote the set of eigen- 
values of matrix P n . Since <£>„ includes eigenvalues —1, 
1, 0, and those derived from 0, $„ can be classified 
into n + 1 non-crossing subsets, denoted by $n,ij 
$n,2i - i $n,n-i, and respectively. $ n ,o contains 

eigenvalues —1 and 1 with single degeneracy; & n ,i con- 
tains all eigenvalues 0; $ n ,2 contains all child eigenvalues 
directly generated from $„ !3 contains all child 

eigenvalues directly generated from $„-i,2; and so on. 

Let Xt = (? ~ 1)( 2g-i +1 " +l? - Then > $« can be explicitly 
expressed as follows: 



where 



and 



$n = (J $n,i , 



$n,0 = {(COSO, 1) , (COS7T, 1)} 



$n,i = (J j ( C0S 
fe=l ^ ^ 



2fc- 1 



-7T, X". 



(C4) 



(C5) 



(C6) 



for 1 < i < n. 

For example, for n — 3, we have 



where 



$3 = $3,0 (J $3,1 |J $3,2 (J $3,3 
$3,0 = {(cosO, 1),(C0S7T, 1)} , 



$ 3,1 = ^(C0S-, — 



$3,2 = 

and 

$3,3 = 



tt (g-l)(2g) 2 +g 
4 ' 2g - 1 



3tt (g-l)(2g) 2 +g 
4 ' 2g - 1 



5tt \ / Tit 



3tt 

COS -,qj , I COS— ,9 j . J cos— .g j . [ COS— .g 



Appendix D: Number of spanning trees in H(q,n) 

Let N st (q,n) denote the number of spanning trees in 
H(q,n). According to the well-known result o 48 ' 60 , the 
number of spanning trees, N st (q,n)^ can be determined 
by N n — 1 nonzero eigenvalues of the normalized Lapla- 
cian matrix for H (q, n), as 



J^dj(n) J^CTj(n) 

N st (q,n) = l=1 - 1=2 . (Dl) 



E^(«) 



1=1 
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Let 6„, A„ and A„ denote Yh=i di(n), JltTi di{n) and 



rii=2 °i ( n ) ' respectively. Then, we have 



and 



N Bt (q,n) 



N st (q,n-1) 



A n -1 X A„_! 



(D2) 



(D3) 



We next determine the recursive relation between 
N s t(q,n) an< i N st (q,n — 1). We first determine the re- 
cursive relations for 0„, A„ and A n . It is evident that 

e„ = J2 = 2£; « = 2 ( 2 <?)" - 2 g e n _!. (D4) 



Moreover, according to the obtained results in the main 
text, it is easy to derive the following recursive relations: 

A„ = ? JV »- 1 A n _ 1 x 2 9B »- 1 = g^- 1 x 2<^"~ 1 A„-i 

(D5) 

and 



A, 



JV«_i-2 



A, 



(D6) 



Combining the above-obtained relations, we obtain the 
recursive relation for N st (q,n) as 

N st (q, n) = q^-^^-^-^N^n- 1) . 

(D7) 

Using N st (q, 0) = 1, we can now solve Eq. (|D7[) to obtain 
the exact expression for the number of spanning trees in 
H(q, n) as 



N st (q,n)=2~ 



2,-^<;-2g 2 +2„g 2 -2" + 1 g " + 1 + 2"+ 1 t ," + 2 



-9 + rag -2ng 2 +2" g " + 1 | 

, «„-n 2 h ™ 



(D8) 



which is in complete agreement with that obtained in— 
yet from a different approach, implying that the tech- 
nique and process of our computation on the eigenvalues 
and their degeneracies for H(q,n) are correct. 
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